Astronomy & Astrophysics manuscript no. pks2005'mwl 


© ESO 2009 


December 24, 2009 





PKS 2005-489 at VHE: Four years of monitoring with HESS and 
simultaneous multi-wavelength observations 

HESS Collaboration, F. Acero 15 , F. Aharonian 1 ' 13 , A.G. Akhperjanian 2 , G. Anton 16 , U. Barres de Almeida 8 *, 

A. R. Bazer-Bachi 3 , Y. Becherini 12 , B. Behera 14 , W. Benbow 1 **, K. Bernlohr 1,5 , A. Bochow 1 , C. Boisson 6 , 
J. Bolmont 19 , V. Borrel 3 , J. Brucker 16 , F. Brun 19 , P. Brun 7 , R. Buhler 1 , T. Bulik 29 , 1. Busching 9 , T. Boutelier 17 , 

P.M. Chadwick 8 , A. Charbonnier 19 , R.C.G. Chaves 1 , A. Cheesebrough 8 , L.-M. Chounet 10 , A.C. Clapson 1 , 
G. Coignet 11 , L. Costamante 1 ' 30 M. Dalton 5 , M.K. Daniel 8 , 1.D. Davids 22 ' 9 , B. Degrange 10 , C. Deil 1 , 
H.J. Dickinson 8 , A. Djannati-Atai' 12 , W. Domainko 1 , L.O'C. Drury 13 , F. Dubois 11 , G. Dubus 17 , J. Dyks 24 , M. Dyrda 28 , 
K. Egberts 1 , P. Eger 16 , P. Espigat 12 , L. Fallon 13 , C. Farnier 15 , S. Fegan 10 , F. Feinstein 15 , A. Fiasson 11 , A. Forster 1 , 
G. Fontaine 10 , M. FuBling 5 , S. Gabici 13 , Y.A. Gallant 15 , L. Gerard 12 , D. Gerbig 21 , B. Giebels 10 , J.F Glicenstein 7 , 

B. Gluck 16 , P. Goret 7 , D. Goring 16 , M. Hauser 14 , S. Heinz 16 , G. Heinzelmann 4 , G. Henri 17 , G. Hermann 1 , 
J.A. Hinton 25 , A. Hoffmann 18 , W. Hofmann 1 , P. Hofverberg 1 , M. Holleran 9 , S. Hoppe 1 , D. Horns 4 , 

A. Jacholkowska 19 , O.C. de Jager 9 , C. Jahn 16 , 1. Jung 16 , K. Katarzyriski 27 , U. Katz 16 , S. Kaufmann 14 , 
M. Kerschhaggl 5 , D. Khangulyan 1 , B. Khelifi 10 , D. Keogh 8 , D. Klochkov 18 , W. Kluzniak 24 , T. Kneiske 4 , Nu. Komin 7 , 
K. Kosack 1 , R. Kossakowski 11 , G. Lamanna 11 , J. -P. Lenain 6 , T. Lohse 5 , V. Marandon 12 , O. Martineau-Huynh 19 , 

A. Marcowith 15 , J. Masbou 11 , D. Maurin 19 , T.J.L. McComb 8 , M.C. Medina 6 , J. Mehault 15 , R. Moderski 24 , 
E. Moulin 7 , M. Naumann-Godo 10 , M. de Naurois 19 , D. Nedbal 20 , D. Nekrassov 1 , B. Nicholas 26 , J. Niemiec 28 , 
S.J. Nolan 8 , S. Ohm 1 , J-F. Olive 3 , E. de Ona Wilhelmi 1 , K.J. Orford 8 , M. Ostrowski 23 , M. Panter 1 , M. Paz Arribas 5 , 
G. Pedaletti 14 , G. Pelletier 17 , P-O. Petrucci 17 , S. Pita 12 , G. Puhlhofer 18 - 14 , M. Punch 12 , A. Quirrenbach 14 , 
', B.C. Raubenheimer 9 , M. Raue 1 ' 30 , S.M. Rayner 8 , M. Renaud 12 ' 1 , F. Rieger 1 - 30 , J. Ripken 4 , L. Rob 20 , S. Rosier-Lees 11 , 
§ 1 G. Rowell 26 , B. Rudak 24 , C.B. Rulten 8 , J. Ruppel 21 , V. Sahakian 2 , A. Santangelo 18 , R. Schlickeiser 21 , FM. Schock 16 , 
U. Schwanke 5 , S. Schwarzburg 18 , S. Schwemmer 14 , A. Shalchi 21 , M. Sikora 24 , J.L. Skilton 25 , H. Sol 6 , L. Stawarz 23 , 
R. Steenkamp 22 , C. Stegmann 16 , F. Stinzing 16 , G. Superina 10 , A. Szostek 23 ' 17 , PH. Tarn 14 , J.-P. Tavernet 19 , 
R. Terrier 12 , O. Tibolla 1 , M. Tluczykont 4 , C. van Eldik 1 , G. Vasileiadis 15 , C. Venter 9 , L. Venter 6 , J.P Vialle 11 , 
P. Vincent 19 , M. Vivier 7 , H.J. Volk 1 , F. Volpe 1 , S.J. Wagner 14 , M. Ward 8 , A.A. Zdziarski 24 , and A. Zech 6 



OS 

o 
o 

(N 

o 

Q 



(N 

o 

u 

43 

Oh 



> 

o> 

O 

(N 



Os 
O 



(Affiliations can be found after the references) 
Received; accepted 

ABSTRACT 



Aims. Our aim is to study the very high energy (VHE; E > 100 GeV) 7-ray emission from BL Lac objects and the evolution in time of their 
broad-band spectral energy distribution (SED). 

Methods. VHE observations of the high-frequency peaked BLLac object PKS 2005—489 were made with the High Energy Stereoscopic System 
(HESS) from 2004 through 2007. Three simultaneous multi-wavelength campaigns at lower energies were performed during the HESS data taking, 
■ consisting of several individual pointings with the XMM-Newton and RXTE satellites. 
Results. A strong VHE signal, ~17er total, from PKS 2005—489 was detected during the four years of HESS observations (90.3 hrs live time). 
The integral flux above the average analysis threshold of 400 GeV is ~3% of the flux observed from the Crab Nebula and varies weakly on time 
scales from days to years. The average VHE spectrum measured from ~300 GeV to ~5 TeV is characterized by a power law with a photon index, 
T = 3.20±0.16 s tat ± 0.10 S yat ■ At X-ray energies the flux is observed to vary by more than an order of magnitude between 2004 and 2005. Strong 
changes in the X-ray spectrum (ATx ~ 0.7) are also observed, which appear to be mirrored in the VHE band. 

Conclusions. The SED of PKS 2005—489, constructed for the first time with contemporaneous data on both humps, shows significant evolution. 
The large flux variations in the X-ray band, coupled with weak or no variations in the VHE band and a similar spectral behavior, suggest the 
emergence of a new, separate, harder emission component in September 2005. 

Key words. Galaxies: active - BL Lacertae objects: Individual: PKS 2005—489 - Gamma rays: observations 



Send offprint requests to: luigi.costamante@stanford.edu or wben- 
bow @ cfa.harvard.edu 
* supported by CAPES Foundation, Ministry of Education of Brazil 
** now at Harvard-Smithsonian Center for Astrophysics, Cambridge, 
USA 

now at W.W. Hansen Experimental Physics Laboratory & Kavli 
Institute for Particle Astrophysics and Cosmology, Stanford University, 
Stanford, USA 



1. Introduction 

PKS 2005-489 is one of the brightest BL Lac objects, at all 
wavelengths, in the Southern Hemisphere. It was initially dis- 
cove red as a strong radio source in the Parkes 2.7 GHz sur- 
vev dWall et al.l 1 1 975l) and later identified as a BL Lac object 
(IWallet alJll98aT It belongs to the complete 1-Jy radio cata- 
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log dStickel et al.ll 19911) . and is o ne of the few extraga lactic ob- 
jects detected in the EUV band dMarshall et all 1 19951) . Its red- 
shift, z — 0.071 dFalomo et al.ll 19871) . is determined from weak, 
narrow emission lines observed during a low optical state. 

PKS 2005—489 is classified as a High-frequ ency peaked 
BL Lac object (HBL; iGiommi & Padovani 1994), because of 
its high X-ray-to-radio flux ratio dSambruna et alj[l995b and be- 
cause its broad-band spectral energy distribution (SED) peaks 
in the UV-soft X-ray band. As is typical of HBLs, the X-ray 
spectrum is dominated by the synchrotron emission of high- 
energy electrons. The second SED hump is expected to peak 
in the GeV— TeV 7-ray band, and is commonly believed to be 
produced by the same electrons up-scattering via the inverse 
Compton mechanism seed photons of lower energy. 

In the X-ray band, PKS 2005-489 has been studied ex- 
tensively, showing an extreme flux and spectral variability. 
Large flux variations with correlated spectral hardening of the 
generally steep spectrum (photon indejfl T = 2.7 — 3.1) 
were o bserved during five EX OS AT observations (Giomm i et alj 
1990; ISambruna etal] |1994|) . Two ROSAT observations in 
1992 confirmed the EXOS AT results and simila rly show a 
soft spectrum (T ~ 3; Sam bruna et alj 1 1995b . A harder 
spectrum (T = 2.3 from 2 to 10 keV) was observed in 
September 1996, during BeppoS AX observations of a brighter 
X-ray state dPadovani et alj|200lh . In October-November 1998 
PKS 2005— 489 underwent a period of exceptional activity, 
with s everal strong X-ray flares. RXTE monitoring observa- 
tions (Perl manet all 1 19991) were performed during the entire 
epoch, and the 2—10 keV flux reached ^3xl0~ 10 ergs cm -2 
s _1 , approximately 30 times higher than the ROSAT values. 
These RXTE observations yielded a detection up to 40 keV and 
showed variations in the photon index between r = 2.3 and 2.8. 
An X-ray flare alert also triggered BeppoSAX observations on 
November 1—2, 1998. The X-ray spectrum was measured be- 
tween 0. 1 to 200 ke V and was characterized by a curved shape 
and harder photon indices, in both the s oft (Ti = 2.0) and har d 
(T 2 = 2.2) X-ray bands (< and >2 keV; lTagliaferri et alj200ll) . 
More recently, observations with the Swift satellite have gener- 
ally sh own PKS 2005—489 in a low-flux, steep-spectrum state 
(r ~ 3. lMassaroet"ai1l2008l) . 

Strong correlations between X-ray an d 7-ray emission 
have been ob s erved in many HBL (e . g. iPian et alj ll 998; 
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. .19991; iKrawczvnski et all l2004t iFossati et alj 
Aharonian et al. 1 120091) . and are typically expected in 
a synchro tron-Compton scenario. Therefore the very bright 
flux and large flux/spectral variability at X-ray energies make 
PKS 2005—489 one of the most promising targets for observing 
a similar behavior in the 7-ray domain. 

In the VHE (E > 100 GeV) band, PKS 2005-489 
was detected neither during observations made between 1993 
and 2000 by eithe r the CANGAROO or Durham groups 



(Rob erts etal J I1998L Il999t iNishiiimal 120021: IChadwick et all 
2000), nor by HESS with a partial array during its commission- 
ing in 2003. In 2004, HESS discovered VHE 7-ray emission 
from PKS 2005—489 with a signi ficance of 6.7a, at a fl ux of 
a few percent of the Crab Nebula (Aharoni an et al. 2005). The 
measured spectrum was soft (T = 4.0). In the MeV-GeV band, 
PKS 2005-489 is one of the few HBL detected by EGRET. 
Howe ver, the observ ed significance is margin al: 3.7a above 100 
MeV dLin etal J 19961) and 4. 1 a above 1 GeV dLamb & Macombl 
119971) . It was instead detected by Fermi with high significance 



2600 bi 
2400 

2200 
2000 

1800 



• On-source 
Off-source 
PSF 




0.05 0.1 0.15 

e 2 [ deg 2 ] 

Fig. 1. Distribution of 8 2 for on-source events (points with statis- 
tical errors) and normalized off-source events (shaded) from ob- 
servations of PKS 2005—489. The curve represents the 9 2 distri- 
bution expected from simulations of a 7-ray point source (pho- 
ton index T = 3.20) at a zenith angle of 40°. The dashed line 
represents the cut on 8 2 applied to the data. 




The power-law spectrum is described as N(E) = No E 
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Fig. 2. Integral flux, I(>400 GeV), measured by HESS from 
PKS 2005— 489 during each dark period of observations (i.e. 
moonless night time within a month). Only the statistical errors 
are shown. The horizontal line represents the average flux for all 
the HESS observations. The four horizontal line segments repre- 
sent the average annual flux observed during the corresponding 
years (see Table[T|i. 



(>1 0q), during the firs t three months of operation (2008, Aug- 
Oct. lAbdo et aDl2009h . 

Because of the high potential for strong 7-ray activity, 
PKS 2005-489 has been monitored at VHE by HESS every year 
since its detection in 2004. Several campaigns of coordinated ob- 
servations with the X-ray satellites XMM-Newton and RXTE 
were also performed. These simultaneous observations, sam- 
pling the same particle distribution through two different emis- 
sion processes, represent a powerful diagnostic tool for probing 
the co nditions of the inner blazar jet, especially during flaring 
events (Coppi & Aharonianl 1999b . 

In this article the results of all HESS observations taken from 
2004 through 2007 are presented, together with the results of 
the multi-wavelength observations. These campaigns character- 
ize the SED of PKS 2005-489 during different states and, for 
the first time, sample both humps of the SED simultaneously. A 
re-analysis of the 2004 HESS data is also provided, which ben- 
efits from an improved calibration of the absolute energy scale 
with respect to the previously published result. 
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Fig. 3. Integral flux, I(>400 GeV), measured by HESS from PKS 2005-489 during each night of observations. The individual plots 
represent each of the four years (2004—2007) of data taking. Only the statistical errors are shown. The horizontal lines represent the 
average annual flux observed during the respective year. The vertical lines at MJD 53282, 53640, and 53642 represent the nights of 
XMM-JVewtoB observations. The epoch of the RXTE observations is between the vertical lines at MJD 53609 and 53623. 
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Table 1. Results from long-term HESS observations of PKS 2005-489. 



Dark 


MJD 


MJD 


Time 


On 


Off 


Q 


Excess 


Sig 


I(>400 GeV^ 


Crab^ 


2 /"KTriTTi 

X (NDF_) 


P(x 2 ) c 


Period 


First 


Last 


[hrs] 










[<t] 


[10 cm s ] 


% 






06/2004 

\J\J/ X-V/ V/ 1 


53171 


53185 


8.1 


678 


5877 


0.0919 


138 


5.4 


2 64 ± 50 


3.0 


7.2 (10) 


0.71 


07/2004 


53199 


53205 


1.5 


105 


985 


0.0909 


15 


1.5 


1.06 ± 1.18 


1.2 


1 4 CD 


0.75 


09/2004 


53255 


53268 


5.6 


342 


3171 


0.0916 


52 


2.8 


1.89 ± 0.78 


2.1 


2 1 (51 


0.84 


10/2004 

1 V.// Z-Y / V./ 1^ 


53285 


53292 


9.0 


569 


5065 


0.0916 


105 


4.5 


2.58 ± 0.62 


2.9 


1 5 f4) 


0.83 


07/2005 


53582 


53595 


9.4 


573 


4504 


0.0908 


164 


7.3 


3.56 ± 0.67 


4.0 


9 2 (9) 


0.42 


08/2005 


53609 


53618 


5.1 


286 


2159 


0.0989 


73 


4.5 


2.86 ± 0.93 


3.2 


4 6 (4) 


0.33 


09/2005 


53639 


53646 


18.1 


1072 


9125 


0.0928 


225 


7.1 


2.93 ± 0.46 


3.3 


13.1 (6) 


0.041 


06/2006 


53908 


53909 


1.3 


127 


824 


0.0938 


50 


4.9 


4.72 ± 1.32 


5.3 


0.0 (1) 


0.97 


a "7 nn/i/c 

07/2006 


53938 


53940 


4.4 


397 


2927 


0.0901 


133 


7.3 


A no LA ^7A 

4.33 ± 0.70 


4.8 


1.3 (2) 


0.52 


08/2006 


53967 


53977 


8.4 


500 


4261 


0.0913 


1 1 1 


5. 1 


2.74 ± 0.57 


3.1 


15.1 (7) 


0.035 


09/2006 


53995 


54002 


7.4 


405 


41 16 


0.0920 


26 


1.3 


1.31 ± 0.58 


1.5 


6.0 (5) 


0.30 


06/2007 


54264 


54270 


5.3 


333 


3006 


0.0924 


55 


3.1 


1.02 ±0.60 


1.1 


1.9 (6) 


0.93 


07/2007 


54291 


54304 


4.4 


309 


2412 


0.0964 


76 


4.5 


2.51 ±0.70 


2.8 


15.2 (8) 


0.056 


08/2007 


54321 


54329 


1.8 


93 


849 


0.0977 


10 


1.0 


0.52 ± 0.92 


0.6 


2.9 (3) 


0.41 


09/2007 


54345 


54345 


0.5 


11 


114 


0.1000 





-0.1 


< 6.60^ 


<7.4 






2004 


53171 


53292 


24.2 


1694 


15098 


0.0917 


310 


7.7 


2.37 ±0.33 


2.6 


2.0 (3) 


0.57 


2005 


53582 


53646 


32.6 


1931 


15785 


0.0930 


462 


11.0 


3.09 ± 0.35 


3.5 


0.6 (2) 


0.73 


2006 


53908 


54002 


21.5 


1429 


12128 


0.0914 


320 


8.8 


2.84 ±0.34 


3.2 


13.5 (3) 


0.0037 


2007 


54264 


54345 


12.0 


746 


6381 


0.0947 


141 


5.3 


1.50 ±0.40 


1.7 


3.8 (3) 


0.28 


Total 


53171 


54345 


90.3 


5800 


49392 


0.0924 


1233 


16.7 


2.57 ±0.18 


2.9 


10.2 (3) 


0.016 



" The integral flux above 400 GeV is calculated using the excess of events above an energy threshold, which differs slightly from the observed 
excess. The quoted error is statistical only and the 20% systematic error on the observed flux is not shown. 
* The percentage is calculated relative to the HESS Crab Nebula flux above 400 GeV ( Aharonian et al . 2006a). 

1 The x 2 , degrees of freedom (NDF), and corresponding x 2 probability P(x 2 ) are given for fits of a constant to I(>400 GeV) binned nightly 
within a dark period, or monthly within a year, or yearly within the total. 
d The upper limit is calculated at a 99.9% confidence level (Feldman & Cousins 1998). 



2. HESS observations and analysis technique 



PKS 2 005-489 was observed with the HESS array dHintonl 
120041) for a total of 158.0 hours (352 runs of -28 min each) 
from 2004 through 2007. During these observations the array 
tracked a position offset from the blazar by 0.5° in alternating 
directions to enable both on-source observations and simulta- 
neous estimation of the background induced by charged cos- 
mic rays. A total of 207 runs pass the standard HESS data- 
quality selection, yielding an exposure of 90.3 hrs live time at 



a mean zenith angle Z n 



35°. The results presented here 



were generated using the standard HESS calibration me thods 
dAharonian et alj|2004l) and analysis tools dBenbowll2005l) . with 
the standard cuts event-selection criteria (except for the 2007 
spectrum, see Sect 3.2). On-source data were taken from a cir- 
cular region of radius 9 cut =Q.ll° centered on PKS 2005-489, 
and the background (off-s ource data) was estimated using the 
Reflected - Regio n method dBerge et alj|2007f) . Equation (17) in 
iLi &Mal (|1983) was used to calculate the significance of any 
excess. All VHE integral fluxes reported throughout this arti- 
cle were calculated assuming the time-average photon index of 
T = 3 . 20 determined in Sect. CO] 



The PKS 2005— 489 observations presented here span four 
years (2004-2007) of HESS data taking. During this time the 
optical throughput of the instrument decreased, because of the 
degradation of the reflective surfaces of the mirrors and Winston 
cones, as well as accumulation of dust on the optical elements. 
For the entire data sample, the optical efficiency has decreased 
by an average of 28% compared to a newly commissioned in- 
strument, with its mirrors installed in Oct. 2001, Dec. 2002, 
June 2003, and August 2003 on CT3, CT2, CT4, and CT1, re- 



spectively. To minimize the effects of long-term variation in the 
optical efficiency of the HESS array, the estimated energy of 
each event was corrected using the ratio of efficiencies deter- 
mined on a run-wise bas is from simulated and observed muons 
dAharonian et alj | 20 06a). After accounting for the decreasing 
optical throughput of the HESS array, the average energy thresh- 
old of the analysis at Z mcan is 400 GeV. 



3. HESS results 

PKS 2005— 489 was clearly detected in each of the four years 
(2004-2007) that it was observed by HESS. A total of 1233 
excess events, corresponding to a statistical significance of 16.7 
standard deviations (er), was detected from the direction of the 
blazar. The results of the HESS observations are given in TableQ] 
which shows the dead time corrected observation time, the num- 
ber of on and off-source events, the on/off normalization (a), the 
excess and the corresponding statistical significance, for various 
temporal breakdowns of the HESS data sample. 

Figure[T]shows the on-source and normalized off-source dis- 
tributions of the square of the angular difference between the 
reconstructed shower position and the source position (8 2 ) for 
all observations. As can be seen in the figure, the distribution of 
the excess (i.e. the observed signal) is very similar to what is ex- 
pected from a simulated point-source of VHE 7-rays at compa- 
rable zenith angles. The off-source distribution is approximately 
flat in 6 2 , as expected. 

The map of excess counts is well-fit by a two-dimensional 
Gaussian with a centroid located at a J2 ooo = 20 h 09 m 27.0 s ± 
1.5* tat ± 1.3* yst and 5 J2000 = -48°49'52" ± 16^ tat ±20^ yst . As 
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Table 2. Best \ 2 fits of a power-law to the various spectra of PKS 2005—489 measured by HESS. 



Epoch 


E lower 


Euppcr 






r 










7400 






2 

X 


NDF 


P(x 2 ) 




[TeV] 


[TeV] 












[10- 


-ii 


cm -2 s 


-1 


TeV" 1 ] 








2004 (AH05£ 


0.2 


2.7 


3.98 


± 


0.38 s tat 


± 


0.10 syst 


0.74 


± 


0.27 s tat 


± 


0.15syst 


5.6 


7 


0.59 


2004^, 


0.30 


4.0 


3.65 


± 


0.39 stat 


± 


0.10 syst 


1.62 


± 


0.33 s tat 


± 


0.32gyst 


6.7 


7 


0.46 


200~ 


0.32 


5.1 


3.09 


± 


0.22 s tat 


± 


0.10 syst 


1.46 


± 


0.28 s tat 


± 


0.29 syst 


13.3 


6 


0.038 


2006 


0.35 


11 


2.86 


± 


0.20 a tat 


± 


0.10 syst 


1.28 


± 


0.20 at at 


± 


0.26 sy st 


3.7 


7 


0.53 


20071 


0.24 


0.75 


3.48 


± 


0.68 stat 


± 


0.10 Syst 


1.00 


± 


0.22 s tat 


± 


0.20 sy st 


0.6 


2 


0.74 


Total 


0.3 


5.3 


3.20 


± 


0.16 sta t 


± 


0.10 sys t 


1.37 


± 


0.17 stat 


± 


0.27 sys t 


9.9 


8 


0.27 



" This (AH05) is the previously published HESS result ( Aharonian et al. 2005) for 2004, and was not corrected for long-term changes in the 
optical efficiency of the system. 

* The 2004 entry is the very same data as presented in AH05, but with a correction (see text) applied to account for relative changes in the 
optical throughput of the HESS array. The 2005, 2006, 2007, and total entries also have this correction applied. 
c The 2007 spectrum is generated with the spectrum cuts iAharoniarP et al. 2006b), see Sect 3.2 



expected, the fit location of HESS J2009-488 is consistent with 
the positio n (a,7 2 noo = 20j 1 9 m 25 .4 s , c5 J2 ooo = -48°49'54") of 
the blazar (Johnst on et al.l IT995h . The upper limit (99% confi- 
dence level) on the extension of HESS J2009-488 is 1 .0' . 

3.1. VHE Flux 

The observed integral flux above 400 GeV for the entire data set 
is I(>400GeV) = (2.57±0.18 sta t±0.51 syst ) x 10~ 12 cm" 2 s" 1 . 
This corresponds to —2.9% of the flux above 400 GeV from th e 
Crab Nebula, as determined by HESS (lAharonian et al.l l2006a). 
Figures [2] and [3] show the flux measured for each dark period 
and night, respectively. 

The integral flux I(>400 GeV) observed during various 
epochs is reported in Table [TJ together with the \ 2 an d corre- 
sponding probability, P(x 2 X f° r fits of a constant to the data 
when binned by nights within each dark period, by dark pe- 
riods within a year, and by year within the total observations. 
There are weak indications (P(x 2 ) < 0.05) of variability on an- 
nual time scales, monthly time scales in 2006, and nightly time 
scales in the dark periods of September 2005 and August 2006. 
The combined variations on different time scales lead to an in- 
dication of variability in the overall monthly curve (x 2 — 30.5 
for 13 NDF, with P(x 2 )< 0.004, Fig. |2]). Although no variabil- 
ity is seen inside the other epochs (e.g. monthly time scales in 
2004, 2005, or 2007), variations in amplitude comparable to the 
statistical errors cannot be ruled out. 



3.2. VHE spectra 

The time-average photon spectrum for the entire data set is 
shown in Fig. [4] All points in the energy range 300 GeV to 
~5.3 TeV are significant except the last (1.5<r) at ~4.6 TeV. 
The data are well-fit, % 2 of 9.9 for 8 degrees of freedom, by a 
power law (dN/dE = h 00 (E / 400 GeV)~ r ) with a photon index 
T = 3.20 ± 0.16 sta t ± 0.10 syst . Removing the -4.6 TeV point 
does not significantly affect the fit result. No evidence is found 
for significant features, such as a cut-off or break, in the energy 
spectrum, also considering the upper limits at higher energies. 

The time-average spectra measured during each year of data- 
taking are shown in Fig. [5] The results of the best \ 2 fits of 
a power law to these data are shown in Table [2] All spectra 



10" 
10~ 1 
10" 1 
10" 1 
10" 1 



2004 - 2007 Data 
Overall Fit 



It 



0.3 0.4 0.5 



4 5 6 7 8 910 
Energy [ TeV ] 



Fig. 4. Time-average VHE energy spectrum observed from 
PKS 2005-489. The dashed line represents the best x 2 fit of a 
power law to the observed data (see Table 2), and then extrapo- 
lated to 10 TeV. Only the statistical errors ar e shown and the up- 
per li mits are at the 99% confidence level (Fel dman & CousinsI 
1998). 



are generated with standard cuts dBenbowl 120051) . except for 
2007. Given the short exposure, low flux, steep spectral slope, 
and degradation of the optical efficiency, a sp ectrum for 2007 
could only be generated with the spectrum cuts ( Aha ronian et al.l 
l2006bl) . which lower the energy threshold of the analysis at the 
expense of sensitivity at higher energies. In Table 2, the epoch, 
lower and upper energy bounds, photon index, differential flux 
normalization at 400 GeV (/400X X 2 ' degrees of freedom (NDF), 
and x 2 probability P(x 2 ) for each fit are given. The x 2 probabil- 
ity for a fit of a constant to the annual T values is 0.30. 

3.3. Correction to 2004 HESS data 



The difference between the two positions is 25" ± 28" 



The data previously published (lAharonian et al.l l2005) for HESS 
observations of PKS 2005—489 in 2004 were not corrected for 
a 19% average decrease in the optical efficiency until 2004. In 
addition, the early simulations neglected the small gaps between 
Winston cone apertures, the combined effect being an overall 
correction of 22% compared to the initial analysis. Figure|6]illus- 
trates the effect of correcting the energy of individual events for 
the relative optical efficiency of the system, as described earlier. 
Comparing (see Tabled the spectrum for 2004 determined here 
to the previously published one shows significant differences in 
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Fig. 5. Photon spectra observed by HESS from PKS 2005—489 during each year of data taking. Only the statistical errors are shown. 
Each line represents the best \ 2 fit of a power law to the observed data. 



the flux normalization (Z400) but not in the photon index (T). The 
flux measured in 2004 is three times higher than previously pub- 
lishecQ, because of the steep spectrum of the source. It should be 
noted that, according to the model for the decrease in the optical 
efficiency of the HESS system, the simulated values match the 
actual values measured in 2003, therefore the upper limit pub- 
lished for 2003 HESS observations is unchange d. Extrapolating 
the 2 003 upper limit (99% confidence level, lAharonian et al.l 
I2005I) to above 400 GeV (using r = 4.0 from the original 
publication) yields I(>400 GeV) < (0.65 ± 0.13 syst ) x 10~ 12 



4. Coordinated XMM-Newton/RXTE observations 

PKS 2005-489 was observed three ti mes with X MM- 
Newton between 2004 and 2005 dobs-id 02 05920401, 
0304080301 and 0304080401). The satellite pointings were 
scheduled such that simultaneous HESS observations were 
possible. The first XMM-Newton pointing was performed 
on October 4, 2004, and yielded a net exposure of 11 ks. 
Unfortunately poor weather did not allow HESS observations 
on this night. Good quality HESS data were taken on the 
following nights. Two more pre-planned XMM-Newton obser- 
vations were performed on the nights of September 26 and 28, 
2005, with exposures of 21 and 25 ks, respectively. The goal of 
scheduling such temporally clos^j XMM-Afe wton pointings was 
to sample both the shortest variability time scales and spectral 
variations occurrin g on the typical tim e scale of HBL (one to a 
few days; see e.g. Tani hata et alj |2001). It should also be noted 



The flux was reported above a threshold of 200 GeV. Extrapolating 
the earlier result, using V = 4.0 as reported, yields I(>400 GeV) = 
(0.86 ± 0.12 stat ± 0.17 syst ) x 10~ 12 cm" 2 s" 1 . 
4 The blazar was observed in 2 consecutive orbits. 
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Fig. 6. VHE spectrum measured by HESS from PKS 2005-489 
in 2004 compared to the previously published version 
(Aha ronian et alj|2005h . Only the statistical errors are shown. 
Correcting the 2004 data for decreases in the optical efficiency 
of the HESS array results in a three-times higher integral flux 
above 400 GeV. 



that the simultaneous observations had to occur near the end of 
the HESS observing season for PKS 2005-489, because of the 
narrow overlap between the HESS and XMM-Newton visibility 
windows caused by the constraints on the orientation of the 
solar panels of the satellite. 

During the XMM-Newton observations, all three EPIC CCD 
cameras were used along with the thin filter. The PN detector 
(lStriideretal.1120011) was operated in timing mode, i.e. when the 
data from a predefined area of the CCD chip are collapsed into a 
one-dimensional row to be read out at high speed. This mode al- 
lows the sampling of the shortest possible flux and spectral vari- 
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PKS 2005-489 PKS 2005-489 

MOS2 PN 4 Oct 2004 broken power-law MOS2PNMOS1 26 Sept 2005 broken power-law 




channel energy (keV) channel energy (keV) 

Fig. 7. Best-fit data and folded model, plus residuals, of the XMM-Newton PN (upper) and MOS2/1 (lower) spectra of 
PKS 2005—489, in October 4, 2004 (left panel) and September 26, 2005 (right panel). The spectra are fitted with a concave broken 
power-law model plus Galactic absorption. 



Table 3. Fits to the X-ray spectra" in the different epochs. 



Instrument 


Band 


Date 


Model" 


Ti 


-E'brcak 


r 2 


^0.1-2.0 keV 


r 2-10 keV 


X i (NDF) 




keV 








keV 




ergs cm 


-V 1 






XMM Mos2+PN 


0.15-10 


4/10/04 


pi 


3.04 ±0.02 






2.46 xlO" 11 


1.03 xlO" 


12 


1.110(284) 








b.p.l. 


3.06 ±0.02 


2.5 ±0.5 


2.64 ±0.20 


2.77 xlO -11 


1.19 xlO" 


12 


1.048 (282) 


RXTE PCU2 


3-20 


20/8-19/9/05 


pi 






2.90 ±0.18 




7.54 xlO" 


12 


0.79(16) 


XMM Mos+PN 


0.15-10 


26/9/05 


pi 


2.31 ±0.01 






7.55 xlO -11 


1.93 xlO" 


11 


1.138(860) 








b.p.l. 


2.34 ±0.01 


0.79 ±0.11 


2.30 ±0.01 


7.72 xl0~ n 


1.95 xlO" 


11 


1.122(858) 


XMM Mos+PN 


0.15-10 


28/9/05 


pi 


2.37 ±0.01 






8.01 xl0~ n 


1.77 xlO" 


11 


1.509(888) 








b.p.l. 


2.39 ±0.01 


1.78 ±0.15 


2.31 ±0.02 


8.19 xlO" 11 


1.83 xlO" 


11 


1.372 (886) 



a The errors are reported at a 90% confidence level for 1 parameter (A% 2 — 2.71). 

Fits with fixed Galactic ./Vh — 3.93 X 10 20 cm~ 2 ; power-law (p.l.) and broken power-law (b.p.l.) models. Unabsorbed fluxes. 



ations without potential pile -urj^ problems during bright states. 
The two MOS (ITurner et al.ll200ll) instruments were used in dif- 
ferent configurations. MOS1 took data in timing mode in 2004 
and in small-window mode during both 2005 observations. The 
MOS2 camera was operated in small-window mode in 2004 and 
in large-window mode in 2005. Simultaneous observations were 
also performed with the Optical Monitor (OM) onboard XMM- 
Newton, using all photometric filters sequentially. For these data 
the central window was read in fast mode to enable temporal 
studies within the data. The exposure for each filter varied be- 
tween 1800 and 4400 s, ultimately constrained by the overall 
duration of the pointing. 

In August 2005, target-of-opportunity observations of 
PKS 2005—489 with RXTE were triggered based on apparent 
enhanced VHE activity in preliminary analysis of uncalibrated 
HESS data. The satellite pointings were scheduled such that fur- 
ther HESS data could be taken simultaneously. Unfortunately, 
only a limited subset of the HESS data in this epoch passes stan- 
dard quality selection criteria because of poor weather condi- 
tions in Namibia. 



4. 1. XMM-Newton data analysis 

The EPIC data were processed and analyzed with SAS 
v7 . 1 . 0, using the calibration files as of July 2008 and stan- 
dard screening criteria^. A period of ~30 minutes at the end 
of Obs-ID 0304080401 was excluded from the analysis because 
of high background rates. The MOS1 and MOS2 imaging ob- 
servations are slightly to moderately affected by pile-up. These 
effects were explored using epatplot (which computes the 
fractions of single and double pixel events showing if they dif- 
fer from normal values), and are sufficiently suppressed by the 
use of only single-pixel (PATTERN=0) events in the analysis. 
Background events for MOS2 are extracted from an annulus 
around the source with inner radius 125" and outer radius 160". 
For MOS 1 the background information was taken from the same 
source region in blank-sky observations, since there is no region 
free of source photons on the same CCD. A check was also per- 
formed using background regions from the outer CCDs, as sug- 
gested in the SAS Analysis Guide. For the PN, the spectra were 
extracted from a rectangular box 20 pixels wide, centered on the 
source strip and extended along the CCD (in raw pixel coordi- 
nates, 28<RAWX<48). The background events were extracted 
from rows 2<RAWX<18, avoiding the noisier strip on the CCD 
border. 



5 That is, when more than one X-ray photon arrives in one camera 
pixel or adjacent pixels before the CCD is read out. 



see SAS User Guide and CAL-TN-0 18. 
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Table 4. Optical Monitor parameters and source fluxes for the 
three XMM-/Vevvfo« observations. 



Filter 


K 




F* onv 




-^4/10 


-^26/9 


-^28/9 


V 


5430 


1.016 


2.49 xlCT 


16 


11.84 


14.51 


14.42 


B 


4500 


1.293 


1.29 xlCT 


16 


9.64 


12.08 


12.09 


U 


3440 


1.634 


1.94 Xl0~ 


16 


7.48 


9.79 


9.78 


UVW1 


2910 


1.870 


4.76 xl0~ 


16 


6.12 


8.27 


8.14 


UVM2 


2310 


2.827 


2.20 xl0~ 


15 


5.57 


7.74 


7.71 


UVW2 


2120 


3.167 


5.71 x!0~ 


15 






7.08 



The effective wavelength of the filter in A. 

The ratio of A\ e , i.e. the extinction in magnitudes at the effective wavelength of 
th e filter A e , to Ay_ (5500A) obtained from the i nterste llar reddening curve given 
bv lCardelli etafl UM% and updated bv lO'Donnelll fl994l) . 

The conversion factor from rate (Counts s - ) to flux (erg cm - 2 s — X A -1 ) from the 
OM in-orbit calibration (see SAS watch-out page). 

The flux of PKS 2005—489, collected for galactic extinction, measured in the three 
epochs (Oct. 4, 2005; Sept. 26 and 28, 2005). The units are mjy. 



The spectral analysis was performed with XSPEC 
vll.3.2ag, testing different binning schemes with at 
least 50 counts in each new bin. For the 2004 observation, the 
spectrum from the MOS1 instrument is not included since the 
timing mode data are significantly noisier than, but consistent 
with, the other detectors. The MOS and PN spectra are fit 
together, with a free constant to allow for different MOS/PN 
normalizations. The PN flux is adopted as the reference value, 
but the MOS fluxes are typically within a few percent of 
the PN fluxes. The X-ray spectra are each fit with single 
and broken power-la w models, with i nterstellar absorption 
modelled by TBabs dWilms et al.l 120001). Th i s mod el is used 
with the solar abundances of Wil ms et al. | (I2000h and the 
cross-sections b y Verner et al. (for a discussion, see 

Baumgartner & Mushotzkv 2006). The results of the fits are 
shown in Fig.|7]and Table [3] (details in Sect. 5.2). 

The data from the OM were processed with the tasks 
omichain for the photometry and omf chain for the timing 
analysis. Since no variations are found in the light curves, the 
whole exposure is used to derive the flux measurements. The 
photometry was processed interactively with omsource for ev- 
ery filter, to assure the use of point-source analysis procedures. 
A standard aperture of 6 pixels on the 2x2 binned images (12 for 
the unbinned images of the 2004 October data set) is used for all 
filters. This corresponds to an aperture of 6" for the optical filters 
(V, B, U). The counts for the UV filters are extrapolated by the 
software, using the UV point-spread function, to an aperture of 
17".5. These are the two apertures for which the OM count-rate- 
to-flux conversion is calibrated. The OM images are affected by 
stray light from a bright star in the field of view, which increases 
the background near one side of the source region. Therefore 
background events are taken from both an annulus around the 
source and from two different circular regions at the same dis- 
tance of the source from the main stray-light reflex. Tests per- 
formed with different background sizes and locations show that 
the photometry does not change by more than 1 %, well below 
the systematic uncertainties of the flux conversion (estimated at 
-10%, see XMM-SOC-CAL-TN-0019). 

The source fluxes were obtained from the count rates using 
the OM conversion factors for white dwarf fl, adding the 10% 
systematic error in quadrature. The fluxes were d e-reddened for 
Galact ic absorption using the extinction curve b v lCardelli et all 
(1989) with the updates by lO'Donnelj (1 1994 . and assuming 
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Fig. 8. Simultaneous X-ray and VHE data taken by RXTE, 
XMM-Newton and HESS in August-October 2005. Upper 
panel: the X-ray photon index. The XMM-Newton values were 
obtained from fits above 2 keV, to be consistent with the RXTE 
energy range. Lower panel: the observed X-ray and VHE flux. 
The X-ray data (filled circles) are energy fluxes integrated in the 
2-10 keV band. The HESS data (open squares) are integral pho- 
ton fluxes above 400 GeV The vertical scales on the left and 
right are for the X-ray and HESS data, respectively. All errors 
are shown at the ler confidence level. 



Rv[= Ay jE(B — V)] = 3.1. This is the average value for the 
Galactic diffuse ISM. For the line of sight of PKS 2005-489, a 
value of A B = 0.241 is used (from NED; ISchlegel et al.lll998l) . 
corresponding to A V (55oo) = 0.182. The conversion factors, ex- 
tinction ratios and the resulting source fluxes are reported for 
each of the filters in Table [4] 

PKS 2005-489 is hosted by a giant elliptical galaxy of total 
R magnitude 14.5 and half-light . radius 5.7" (from HST snapshot 
observations. IScarpa et all2000l) . For the SEDs (Figs. 191 andlTUl, 
the OM fluxes were corrected for the contribution of the host 
galaxy. The wavelength-dependent correctio n was determine d 
using a template SED for elliptical galaxies (ISilva et al.| [l998). 
rescaled to the host-galaxy flux in the R band, and accounts for 
the given apertures. The OM fluxes are always dominated by the 
non-thermal emission, and the small contribution of the host- 
galaxy is only noticeable in the V and B filters. 



4.2. RXTE data analysis 



http://xmm2.esac. esa.int/sas/7. 1 .0/watchout/ 



The RossiXTE/PCA (Jahod a et all Il996h performed 15 snap- 
shot observations of PKS 2005— 489 between August 20 and 
September 19, 2005, yielding a total exposure of 23.4 ks. The 
observations were mostly done with PCU0 and PCU2. The PCA 
STANDARD2 data were reduced and analyzed with the standard 
FTOOLS routines in HEASOFT V6.3.2, using the filtering cri- 
teria recommended by the RXTE Guest Observer Facility. The 
new SAA history file and parameters for the background calcu- 
lation provided in September 2007 were used, which correct a 
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bug in the calculation of the faint background modefl Only the 
top-layer events were processed and only the PCU2 data were 
considered, to ensure a more accurate spectral measurement. 
The average net count rate in the 3-15 keV band is 0.54 ± 0.02 
cts s _1 pcu -1 . The RXTE spectra were extracted and fitted sep- 
arately for each pointing, and summed together to obtain the av- 
erage spectrum. Each spectrum is well-fit by a single power-law 
model. 



5. XMM-Newton/RXTE results 

5. 1. Flux variability 

The 2004 and both 2005 EPIC light curves contain no evidence 
of flux variability within each exposure, and have average net 
count rates in the MOS2 camera of 1.14 ± 0.04, 6.38 ± 0.08, and 
6.15 ± 0.08 counts s _1 , respectively. The probability of constant 
emission is P CO nst > 0.99, 0.92 and 0.99 for the three epochs, 
respectively. Similarly, no flux variability is observed in the OM 
data for any of the observed filters. No significant flux variations 
(Pconst > 0.76) are detected during the RXTE observations (see 
Fig. [8]l as well. However, night-by-night flux variations of a few 
tens of percent cannot be excluded, due to low statistics in the 
RXTE measurements. 

Although no variability is found within any of the XMM- 
Newton exposures or within the complete RXTE sample, strong 
flux variations are found on longer time scales. The 2-10 keV 
flux varies by a factor of ~46 between 2004 and 2005. In the 
UV and optical bands the flux increases by ^40% and ^20%, re- 
spectively, between 2004 and 2005. A smaller (factor of 2.5), but 
significant, change in the 2-10 keV flux also occurs during the 17 
days between the final RXTE and first XMM-Newton pointing 
in September 2005. 

5.2. X-ray spectrum 

Since there is no variability in the X-ray flux or hardness ratio, 
during each XMM-Newton exposure, a single X-ray spectrum 
was extracted for each of the three pointings (see Table [3] and 
Fig. [7]l. The spectra from the individual RXTE exposures do not 
vary significantly (see Fig. [8} and are consistent with the fit to 
the complete RXTE sample reported in Table [3] In each case 
(XMM-Newton and RXTE), the Galactic absorption is fixed to 
A H = 3.93 x 10 20 cm" 2 , as recently determined from the LAB 
survey (Kalberla et al. 2005). It is important to note that there is 
a discrepancy between the L AB survey estimate and the value 
obtained from the HI maps bv lDickev & Lockmanl (fl990) (DL, 
Ah = 5.08 x 10 20 cm -2 ). The LAB survey value is adopted 
because it provides better residuals, is consistent with the re- 
sults when A H is left as free parameter ((3.96 ± 0.17) x 10 20 
cm -2 ), and is close to the value previously adopted for the 
BeppoSAX results on the large flare of 1998 (4.2 x 10 20 cm -2 , 
Tagliaferr i et al.l2001h . Using the DL value significantly changes 
the slope of the source spectrum below 1 keV. Therefore both 
cases are considered for interpreting the fits results, as additional 
systematic uncertainty. 

In all three data sets, the EPIC spectra are better fit by a bro- 
ken power-law model with respect to a single power law, with 
high significance (F-test > 99.9%). Remarkably, the parameters 
reveal an inverted broken power-law spectrum, where the slope 
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Details at http://www.universe.nasa.gov/xrays/programs/rxte/pca/ 
doc/bkg/bkg-2007-saa/ 



Fig. 9. SED of PKS 2005-489 in the frequency range covered 
by the XMM-Newton and RXTE observations (filled circles). 
Archival X-ray spectra are shown along with those measured by 
XMM-Newton on Oct. 4, 2004 (red points), RXTE in Aug-Sept 
2005 (green points) and XMM-Newton on Sept. 26, 2005 (blue 
points). For the latter, two fits are plotted, showing the effects 
of using two different estimates of the Galactic column density 
(A H = 5.08 or 3.93 x 10 20 cm~ 2 ; see text). For better visibility, 
the XMM-Newton spectrum of Sept. 28, 2005 is not plotted as 
it would appear superposed on the Sept 26 data. The flux of the 
host galaxy in the R band is plotted as black square. Historica l 
data are shown in gray (from NED and Tagliaf erri et alj|2~001l) . 
In the X-ray range, they correspond to (from top to bottom ): the 
BeppoSAX data of th e 1998 flare (Tagliaf erri etai] 1200 ll) . the 
BeppoSAX 1996 data dPadovani et alJl2001l) . the S wift spectrum 
in April 2005 (filled bow-tie. iMassaro et al.ll2008l). and the 1992 
ROSAT spectrum (open bow-tie. IComastri et 19971) . 



in the soft X-ray band is steeper than in the hard X-ray band 

(r 2 <r!). 

In the 2004 dataset, the X-ray spectrum is characterized by 
the lowest flux ever recorded from this source (see Fig. [9]) and a 
very steep slope of T ~ 3 up to ^3 keV. At higher energies, there 
is evidence of a moderate hardening trend. The trend is consis- 
tent between the MOS and PN spectra, but is only significant in 
the PN data thanks to better statistics. This hardening can also 
be reproduced by the sum of two power-law models, the second 
having a harder index. The slope of this second power-law func- 
tion is not well constrained. However, assuming a fixed value of 
r = 1.5, as could be expected for the inverse Compton emission 
of low-energy electrons, a good fit (x 2 = 1-07 for 283 NDF) 
is obtained with a flux normalization at 1 keV for the second 
power-law function at 2.8 ± 0.7% of the first power law. In this 
case, the two power-law functions would cross at about 20 keV. 
Using the high Ah has minimal impact on this result, steepening 
the low-energy slope by only +0.1. 

The difference between the EPIC spectra of September 26 
and 28, 2005 is negligible. However, the spectrum from either 
of the 2005 observations is significantly harder and brighter than 
in 2004 (Fig. |9). Either 2005 spectrum is represented well by 
a pure power-law model with T = 2.3 from ~1 to 10 keV. 
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Table 5. Best power-law fits to the HESS spectra c orrected for 
EBL absorption following lFranceschini et all ((2*008). 



Epoch 


■pa 
int 


i^(lTeV) b 
[10~ 12 erg cm" 2 s" 1 ] 


X 2 (NDF) 


2004 
2005 
2006 
2007 


3.1 ±0.4 
2.6 ±0.2 
2.4 ±0.2 
3.1 ±0.7 


1.8 ±0.4 
2.7 ±0.4 

2.9 ±0.4 
1.2 ±0.9 


6.4(7) 
13.2(6) 
4.3 (7) 
0.5(2) 


2005 7-8 c 
2005 9 C 


3.2 ±0.4 
2.5 ±0.3 


2.4 ±0.6 
2.9 ±0.5 


9.2(4) 
4.7(6) 


Total 


2.69 ±0.16 


2.34 ±0.20 


8.9(8) 



a Only the statistical errors are shown. 

b The normalization at 1 TeV, in vF u units. 

c The HESS spectrum in 2005, extracted in two different epochs (see Fig.FJ}: before the 
XMM pointings (observing periods 2005, 7-8) and during the XMM epoch (period 
2005,9). 



At lower energies, each has a slightly steeper slope. Although 
modest (Ar < 0.1), the spectral break is significant (F-test 
> 99.99%) and is not caused by the extension of the MOS data 
at low energies, or by the cross-calibration between the detec- 
tors. The broken power-law model is also statistically required 
for the PN spectrum alone (F-test > 99.99%), as well as for the 
MOS spectra fitted separately. It should be noted however that 
the precise location of the break is affected by systematic more 
than statistical errors. In Table 3, the break energy for the fit with 
the overall minimum \ 2 is reported, but there are other similar 
local minima between 0.7 and 2.2 keV, whose hierarchy in \ 2 
can change according to the calibration of the effective area near 
the instrumental Si and Au edges (at 1.7 and 2.1 keV). The true 
break therefore should be considered more realistically inside 
the range 0.7-2.2 keV. This does not affect the values of the soft 
and hard slopes in a relevant way (they do not change by more 
than 0.04 and 0.02, respectively). 

The amount of the spectral break, however, depends on the 
adopted value of the Galactic column density. Using the DL 
value instead of the LAB survey value, the soft X-ray slope be- 
comes steeper (Ti = 2.6), making the break more pronounced 
(AF ~ 0.3). The two cases are shown in Fig. [9] where the spec- 
tra are plotted together with the other data sets and archival data. 
Allowing instead the column density to be lower than all avail- 
able estimates, a single power-law spectrum is obtained with 
N H = 3.76 ± 0.08 cm" 2 ( X 2 r = 1.226 for 859 NDF), though 
a broken power-law model still seems to provide a better repre- 
sentation (F-test > 95%). 

The 2005 RXTE observations yield a very steep (r = 2.9) 
spectrum from 3-15 keV. This spectrum is almost parallel to 
the 2004 XMM-Newton spectrum but with ~6 times higher 
flux. The RXTE photon index and flux are in g ood agreement 
with t he values measured by Swift in April 2005 dMassaro et al.l 
2008). This suggests that the hardening of the X-ray spectrum of 
PKS 2005 -489 occurred in September 2005, between the RXTE 
and XMM-Newton pointings. 

6. Discussion 

The spectral energy distribution (SED) of PKS 2005-489, as- 
sembled with the data taken in 2004 and 2005, is shown in 
Fig. [10] The VHE 7-ray spectra are corrected for the energy- 
dependent attenuation caused by interactions with the optical- 
IR photons of the diffuse extragalactic background light (EBL). 



The correction w a s per formed adopting the EBL model of 
iFranceschini et al.l (120081) . This model takes into account the 
most recent results on galaxy emission and evolution, and 
agrees with b oth the upper limits on the EBL derived from 
blazar spectra ( Aha ronian et al . 2006c, |2007j) and the lower lim- 
its give n by galaxy counts, at both optical and inf rared wave- 
lengths dMadau & Pozzettill2000t iFazio et al.ir2004l) . The fits to 
the absorption-corrected spectra are reported in Table [5] and 
plotted in Fig.fTOl 

The simultaneous observations do not show evidence of 
strong changes in the location of the SED peaks over the years, 
with respect to the historical values. The hard optical-UV spec- 
trum from the OM photometry and the steep X-ray spectrum 
above 0.1 keV locate the synchrotron peak between the two 
bands, at approximately 0.5 x 10 16 Hz. At VHE, the steep spectra 
constrain the peak of the 7-ray emission to energies <0.2 TeV. 
The low flux in the VHE band does not enable a determination of 
the spectral properties in exactly the same (short) observing win- 
dows of the X-ray observations. However, given the lack of sig- 
nificant VHE variability in 2004 or 2005, the monthly or yearly 
average spectra can be considered a reasonable approximation 
to the VHE spectra during the epoch of the X-ray observations. 

In 2004, PKS 2005 -489 was in a very low state in both 
the X-ray and VHE bands. The spectral slopes are also simi- 
lar (r ~ 3) for the X-ray and VHE emission. The X-ray flux 
around few keV is more than two orders of magnitude less than 
the BeppoSAX flux during the active state of 1998. The slight 
spectral hardening towards higher X-ray energies is consistent 
with the onset of the concavity expected from the emergence of 
the inverse-Compton emission in the SED. This emission, pro- 
duced by low-energy electrons, is usually not visible in the X-ray 
spectra of HBL because of the dominance of the synchrotron 
emission from high-energy electrons. In this data set, its pres- 
ence becomes detectable because of the steepness and very low 
flux of the synchrotron (X-ray) spectrum. A similar behavior has 
recen tly been observed in the HBL PK S 2155-304 during a low 
state (Fos chini et al. 2008; Zhang 2008). This transition between 
steep synchrotron and flat inverse-Compton spectra is instead 
fully visible inside the observed X-ray passband for int ermediate 
BL Lac objects (e.g. ON 231, see lTaghaferri et al.ll2000T) . 

In 2005, the X-ray spectrum changes radically. Swift and 
RXTE measurements made in April and August-September 
2005, respectively, show a large increase in flux but little change 
in the spectral slope compared to 2004. Shortly after the RXTE 
observations, a much harder spectrum (Arx — 0.7) and a 
brighter state than in the previous epochs is found during XMM- 
Newton observations. Most interestingly, the data require a spec- 
tral break with a steeper spectrum towards lower energies. This 
break is more apparent when a higher galactic Nh is used (e.g. 
the DL versus LAB survey values; see Fig. [9), but even when 
considering a single power-law shape with an ad-hoc Nh, the 
extrapolation of the X-ray spectrum down to lower frequencies 
does not match the UV-optical data, which are still rising with 
frequency. A hump with a steeper soft X-ray slope is needed to 
smoothly connect with the optical-UV spectrum, and this feature 
indicates the presence of two particle populations, since it can- 
not be easily reproduced by synchrotron radiation with a single 
electron distribution. 

At VHE, the comparison between the 2004 and 2005 spec- 
tra indicates a hardening trend as well, though with marginal 
significance (~1 sigma; see Table |5]l due to the large errors on 
the HESS measurement. However, if this trend is real, and if 
the VHE spectrum is indeed produced by the same particles 
as are responsible for the X-ray emission (as suggested by the 
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Fig. 10. SED of PKS 2005-489, in differ ent epochs and states. The XMM-Newton a nd RXTE data are plotted as in Fig.[gJ together 
with historical data (shown in grey; from lTagliaferri et alJl200lMMassaro et al . 2008). The recent Fermi-LAT spectrum (0.2<E<10 
GeV) from Aug-Oct 2008 is also plotted, for reference (open bow-tie; lAbdo et aLll2009l) . For clarity, the HESS sp ectra (E>300 
GeV) are plotted as bow-ties, and have been corrected for EBL absorption with the model by lFran ceschini et al. (2008, see Table 5). 
Main panel: the HESS time-averaged spectrum from 2004 (red filled bow-tie) and 2005 (blue hatched bow-tie). Right mini-panel: 
the HESS spectrum in 2005, divided in two subsets: before the XMM pointings (periods 07-08/2005 in Fig.[3j green filled bow-tie) 
and during the XMM epoch (period 09/2005 in Fig.[3j blue hatched bow-tie). Same colors associate optical, X-ray and VHE datasets 
corresponding to same epochs. 



similar spectral slopes in 2004-2005 and by the location of the 
synchrotron pealfl), a similar behavior is expected. Namely, the 
VHE spectral hardening should have occurred after the RXTE 
observations, during the epoch of the XMM-Newton pointings, 
while the 2005 July-August VHE spectrum should be similar to 
the 2004 spectrum (modulo normalization). 

To test this hypothesis, the 2005 HESS data set was di- 
vided into two subsets. The first consists of all VHE data be- 
fore the XMM-Newton pointings (2005 July-August dark pe- 
riods), and the second contains all 2005 HESS data taken dur- 
ing or after the pointings (2005 September dark period). Table|5] 
shows the results of power-law fits to the spectra from these sub- 
sets, after correction for the absorption of VHE photons on the 
EBL. The HESS spectrum in the XMM epoch is indeed harder 
( Ar int = 0.72 ±0.47) than the pre-XMM spectrum, and the sta- 
tistical significance of the spectral change increases (^2er) with 
respect to the comparison between yearly spectra. In addition, 
the pre-XMM spectrum has the same slope as the 2004 HESS 
spectrum. The HESS data therefore are fully consistent with a 



If the VHE electrons were emitting by synchrotron at energies well 
below the X-ray range, they would correspond to the synchrotron peak, 
around 10 16 Hz. Since the SSC scattering at the peak would be in the 
Thomson regime for this frequency, the gamma-ray peak should then 
be inside the VHE band, which is not observed. 



behavior in the VHE band that mirrors the spectral variations 
seen in the X-ray band, as expected if both emissions are pro- 
duced by the same electrons. Interestingly, the harder VHE state 
continues in 2006. Because of the limited exposure in 2007, the 
VHE statistics are too poor to draw any meaningful conclusions. 

The most remarkable feature of the VHE emission, however, 
is the lack of flux variability. Comparing the average states in 
2004 and 2005, the integrated energy flux (e.g. for reference in 
the decade 0.3-3 TeV) does not vary by more than 40% between 
2004 and 2005, and is consistent with a constant flux. In contrast, 
the X-ray flux increased by a factor of ~ 16 in the 2-10 keV band 
(see Table [3). The optical-UV emission, which is close to the 
synchrotron peak and which typically provides the target pho- 
tons for the IC scattering in the VHE band, shows an increase 
as well (^40%). In a homogeneous SSC scenario, if the VHE 
emission is produced by the same electrons emitting in the X- 
ray band (as indicated in this case by the spectral behavior), a 
fresh injection of electrons should make the VHE flux increase 
at least linearly with the X-ray flux. The reason is that the energy 
density of all possible seed photons for the IC scattering has in- 
creased or is constant between the two epochs. For the VHE flux 
to remain constant, a corresponding decrease in the seed photons 
energy density is required. 
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The X-ray and VHE data taken together, therefore, suggest 
that a new jet component is emerging in the SED of PKS 2005- 
489, which is physically separated from the main emitting re- 
gion. The emission at the synchrotron peak of this new com- 
ponent should be lower, remaining hidden below the observed 
SED, while its harder spectrum at high energies emerges in the 
hard-X band. The electrons of this new component would not 
see the energy density of the observed synchrotron peak, but the 
lower energy density of their self -produced synchrotron peak in- 
stead. 

7. Conclusions 

PKS 2005-489 is detected at VHE in each of the four years 
it was being observed by HESS (2004-2007). The 2005-2007 
data clearly confirm the VHE discovery reported by HESS 
dAharonian et al.ll2005l) and quadruple the statistics of the initial 
spectrum measurement. Re-analysis of the previously published 
2004 HESS data, using the improved calibration of the detec- 
tor's energy scale, results in a ^3 times higher flux and a similar 
photon index. 

The measured VHE 7-ray flux is low (~3% Crab) and only 
shows weak variations on time scales ranging from days to 
years. The observed time-averaged (2004-2007) VHE spectrum 
is soft, with a photon index V = 3.20 ± 0.16. Although evi- 
dence of VHE spectral variations is marginal by itself, the VHE 
spectrum seems to track the X-ray slope variations when multi- 
wavelength coverage is available. Observations performed with 
XMM-Newton and RXTE in 2004 and 2005 reveal remarkable 
changes in the X-ray spectrum, but without shifting the location 
of the synchrotron peak with respect to historical observations. 

Interpreting these measurements along with the HESS data 
suggests the emergence of a new jet component in the SED 
that is characterized by a harder electron spectrum. This compo- 
nent must be separate: its particles cannot interact with the syn- 
chrotron photons of the observed SED peak, otherwise higher 
VHE fluxes than observed would be implied. 

PKS 2005—489 is found overall in a very low state, in both 
the X-ray and VHE bands, during the observations presented 
in this article. PKS 2005—489 has historically demonstrated a 
lOOx dynamical range in the X-ray band. Thus, dramatically 
higher VHE fluxes (10 2 — 10 4 x) can be expected in the future, 
unless such an increase in the X-ray flux is counter-balanced by 
a strong (>10x) and simultaneous increase in the blazar's mag- 
netic field. Further monitoring of this object is highly encour- 
aged, as it is one of the few HBL easily detected at VHE during 
low states and has the potential for extreme brightness and vari- 
ability. 

The results presented here confirm the strong diagnostic po- 
tential of coordinated optical-X-ray-VHE observations. Future 
studies can be significantly improved by incorporating data from 
the recently launched Fermi 7-ray satellite. The Fermi data will 
provide information on the lower energy side of the inverse- 
Compton peak, enabling contemporaneous measurements of 
both sides of each blazar hump. 
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